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1 Introduction 

Given a natural Hamiltonian L with n degrees of freedom, satisfying some additional geometric 
conditions, it is shown in [1] how to generate a n + 1 degrees of freedom Hamiltonian H, called 
the extension of L and depending on an integer parameter m G N \ {0}, such that H admits two 
new independent first integrals: H itself and a polynomial in the momenta of degree m. This 
implies that, if L is superintegrable with 2n — 1 independent first integrals, then all the extended 
Hamiltonians H also are superintegrable for any m with the maximal number of 2n + 1 = 2(n + 
1) — 1 first integrals, one of them of arbitrary degree m whose expression is explicitly computed by 
means of a simple iterative process. The extension procedure, summarized in Section 2, is applied 
to the superintegrable systems on (Section 3) and (Section 5) as listed in [5]. These are 
all the superintegrable systems on and ¥? admitting two independent first integrals of degree 
two in the momenta other than the Hamiltonian. It is found that a great part of them admits 
superintegrable extensions and for some of them the extensions are explicitly determined. In 
Section 4 the possible natural Hamiltonians admitting an extension are determined on E" and the 
superintegrable systems of Calogero and Wolfes are considered as examples in E^. In Section 6 
the extension procedure is applied to a class of Hamiltonians on constant curvature manifolds to 
which some generalizations of the Tremblay-Turbiner-Winternitz (TTW) system belong [10]. 
These generalizations of the TTW systems are superintegrable for rational values of certain 
parameters, for which they admit polynomial first integrals of degree greater than two [7, 8]. 
The cases allowing the extensions are determined and the TTW systems are among them, 

*This paper is a contribution to the Special Issue "Superintegrability, Exact Solvability, and Special Functions" . 
The full collection is available at http://www.emis.de/journals/SIGMA/SESSF2012.html 
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obtaining in this way their superintegrable extensions. Other superintegrable generalizations in 
higher dimensions of the TTW system are obtained in [6]. 

2 Extensions of superintegrable systems 

We resume in the following statement the main results proved in [1]. 

Theorem 1. Let Q be a n-dimensional {pseudo-)Riemannian manifold with metric tensor g. 
The natural Hamiltonian L = ^g^-^PiPj + V{q^) on M = T*Q with canonical coordinates {pi,q^) 
admits an extension 

H = ^pl + aiu)L + fiu) (1) 

with a first integral F = U'^{G) where 
U = Pu + ^{u)Xl, 

Xl is the Hamiltonian vector field of L and G{q^), if and only if the following conditions hold: 
i) the functions G and V satisfy 

H(G) + mcgG = 0, mGN\{0}, c G M, (2) 
VF- VG-2m(cy + Lo)G = 0, Lq G M, (3) 

where H(G)jj = ViVjG is the Hessian tensor of G. 
ii) for c = the extended Hamiltonian H is 

H = ]^pl + mA{L + VQ)+mLQA^[u + UQf, (4) 

for c / the extended Hamiltonian H is 
1 2 m(cL + Lo) 

with K, uq, Vo, Wq, a £R, a ^ and 

sin y/KX 



K> 0, 

K 

X, K = 0, 

sinh a/|k|x 

K < 0. 



Snix) 



Dynamically, extended Hamiltonians (1) can be written as 

^Pu - c2t "^1 \^ -h = 0, cL + Lo + v = 0, if c / 0, 

2 S^{cu + uo) 



and 



^pI + mLoA'^u^ - r]-h = 0, mA{L + Vb) + r/ = 0, if c = 0, 
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with H = h, and where constant rj can be understood either as a separation constant, between 
the elements of H depending on {u,pu) and those depending on {q^,Pi), or as a couphng constant 
merging together the Hamiltonian L on T*Q and the Hamiltonian 

^Pu-^( \ ^' or -p^ + mLoAu, c = 0, 

2 bf.{cu + UQ) 2 



depending on {u,pu), to build the extended Hamiltonian H. Several examples are given in [1]. 

dx 



Under the hypotheses of Theorem 1, recalling that Ck(x) = 4^,Sk{x), the polynomials in the 



momenta 

(Pu - A{u + uo)XLr G, c = 0, 
are first integrals of H of degree m. For example, 

UG = Gpu + i{u)Xl{G) = Gpu + l{u){L, G}, 
where { , } are the Poisson brackets. Another way to calculate f7™G is to apply the formula [2] 

C/"^G = PmG + n^XLG, 

with 

[m/2] 



1 — n \ / 



k=0 

[m/2]-l 

E (2fc + i)7'''''^'r^'"'(-Mci> + io)f, m>l, 

A;=0 ^ ^ 

where [•] denotes the integer part, Di = 7 and 

{ C^(cu + Up) ^ _^ ^ 
^^(cu + Mo)' ' (6) 

-A{u + uo), c = 0. 

Remark 1. For c = 0, if Lq = then p„ is a first integral of the extended Hamiltonian and 
the extended potential is merely mAV, with m and A constants. Therefore the extension is 
trivial. Otherwise, a harmonic oscillator term in the variable u, attractive or repulsive, is added 
to the potential mAV . In the case c 7^ 0, we remark that, in the extended Hamiltonian (5), the 
potential V is multiplied by a non-constant factor depending on u. 

Remark 2. The expression (6) of ^{u) is determined in [1] as the general solution of the 
differential equation 

7' + 0(7^ + = 0, (7) 

for real values of 7, u, c and k. However, the general solution of equation (7) in the complex 
case can be considered too. Its expression is the same as (6) if we extend functions S'k and 
to complex values of x and n by using the standard exponential expressions of trigonometric 
functions. After this generalization, the extension procedure characterized by Theorem 1 can 
be applied also to the complex case. 
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In [1] it is proved that, if n > 1, then (2) admits a complete solution G depending on 
a maximal number of parameters (oi), with i = 0, ...,n, iff the sectional curvature of Q is 
constant and equal to mc. Once such a solution G is known, an extension of L is possible iff the 
compatibility condition (3) on V is satisfied. A sharper result on constant curvature manifolds 
is the following 

Proposition 1. On a manifold with constant curvature K the only eigenvalues mc for the 
Hessian equation (2) are either zero or the curvature K. Moreover, if K ^ and mc = then 
the only solutions of (2) are G = const. 

Proof. The equation (2), written in components, becomes 

ViVjG + mcgijG = d^G - V'ljdkG + mcgijG = 0. 

and the integrability conditions that each solution G{q^) must satisfy are given by (see [1] for 
details) 

Rhij^k = 'mcigjhZi- gihZj), V/i, V^t^j, 

where Zk = dkG/G and R'^-j = O^rJ; - djVfi + Tj^rf^ - Tf^rJ^ is the Riemann tensor of the 
metric. For a constant curvature manifold we have Rhuj = K^gjig^i — gughj) and the above 
conditions become 

{K - mc){gjhZi - gihZj) = 0, V/i, ^i^^j- 

By choosing orthogonal coordinates, we see that, since i ^ j, the equations are identically 
satisfied for h ^ i,j. Otherwise, they reduce to {K — mcjgjjZi = 0. Hence, for mc ^ K, the 
only possibility is = for all i (that is, G is a constant). For mc ^ K and mc ^ 0, by (2) we 
get G = 0, thus mc is not an eigenvalue. ■ 

In [1] it is shown that the first integral U'^{G) is functionally independent from L and 
all its possible first integrals Li in T*Q. It is straightforward to see that L and Li are first 
integrals of therefore, if L is a superintegrable Hamiltonian with 2n — 1 first integrals, inclu- 
ding L, then H is superintegrable too with 2n — 1 + 2 = 2(n + 1) — 1 first integrals including H 
itself. It follows that the extension procedure applied to a superintegrable Hamiltonian L, under 
the hypothesis of Theorem 1, always produces a new superintegrable Hamiltonian H. Given 
a superintegrable system of Hamiltonian L = ^g^-'piPj + V{q^) with a configuration manifold Q 
of constant curvature K, its extension to another superintegrable system of Hamiltonian H, 
when possible, can be obtained by applying the following algorithm (see also [1]): 

1. Solve equation (2) on the manifold Q, with c = K/m, to compute the general form of the 
function G(g*, oq, . . • , a„). 

2. Solve equation (3) with the given V for some of the parameters (oj) in G. This is a crucial 
step, because if no solution is found, except for the trivial one G = 0, then the extension 
is not possible. 

3. Determine the extension via Theorem 1. 

4. Compute U"^{G) to obtain the additional first integral. 

In the following sections we analyze several examples of extensions of superintegrable systems. 
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3 Superintegrable extensions of systems 

Since the curvature of is zero, equation (2) admits, by Proposition 1, solutions G ^ only for 
mc = 0. Thus, we are in the case c = of Theorem 1 and the extended Hamiltonian H is in the 
form (4). In [1] the complete solution G of (2) is computed, in standard Cartesian coordinates 
of E^, as 

G = ao + aix + a2y. (8) 
Equation (3) becomes here 

d^yd^G + dyVdyG = 2mLoG, 
whose general solution is 

V = mLo[{x + + (y + 2/°)'] + F{a2X - aw), (9) 

with the constraint oq = aix^ + a2y^, where x^,y^ £ or C and F is any regular function of 
the argument. The extension, after the non restrictive assumptions Vq = uq = and A = m~^, 
becomes 



TT ^ 2 , T , -^0 2 ^ / 2 , 2 , 2\ , , -^0 2 

2 ml ^' m 

In [5] the list of all superintegrable potentials in E^ with three independent quadratic in the 
momenta first integrals is given, up to isometrics and reflections. In that article E^ is assumed 
to be a two-dimensional complex manifold. According to Remark 2 we apply in this case the 
same procedure developed for the real case and allow all functions, variables and parameters 
(except for m G N) to take indifferently real or complex values, in this one and all the following 
sections. By setting z = x + iy, z = x — iy (remark that, despite the notation, if x and y are 
complex coordinates, then z and z are not complex conjugate one of the other), the list is 

El V = ^ + ^+a^{x' + y'), 
E2 y = aix+p+a3(4x2 + y2)^ 

E3 y = a3(x' + y'), 

E4 V = ai{x + iy), 
E5 V = aix, 



E6 V 



ai 



E7 V = = H , , h aszz, 

Vz^ - Vz^-k^iz + Vz^-k^y 

aiz 02 , 
E8 V = -z^ + ^ + a3zz, 

Z'^ Z'^ 

ai x + z 

-C/9 \/ = — = + 022; + "3 — 7^-, 
y z \J z 

( 3_2\ / _ 1_V 

£■10 V = a\z + a2 Iz — -z j + \zz — -z 

Ell V = aiz -\ 7^ + ^, 

V z \J z 



E12 V 



aiz 
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E13 V 



ai 



EU V 

Elb V 

EIQ V 

En V 

El^ V 

E19 V 

E20 V 



h(z), for any function h, 



1 



a\ + 



OL2 



X 



+ 



as 



03 



ai ao 

-T^ + ^+ , - 
y zz ^ z\J zz 



\/ x"^ + y'^ 
aiz 



+ 



a2 



+ 



as 



1 



ai + a2\/ X + + y2 _|_ 



Y''x - \/x2 _|_ y2 



with (aj), A; € C. The equation of compatibihty (3) is equivalent to a Hnear homogeneous 
expression in (oj) with coefficients hnear but not homogeneous in a^. This expression vanishes 
only for some suitable choices of the parameters Oj and a^, depending on the extension parame- 
ters m, Lq. The solution (9) of (3) shows that a non null function G of the form (8) satisfying (3) 
exists only for the following potentials 



V 



G 


particular cases of 


aix + a2y 


El, E3, E7, E8 




El 


aix 


El 




E2 




E2 


aiz 


E7 


aiz 


E8 




ElO 



I 

ii 
Hi 
iv 

V 

vi 
vii 

via 



mLo{x'^ + y"^) 

^ + mLo{x'^ + 2/2) 

^ + mLoix'^ + y"^) 

y2 

aix + mLo(4x2 + y^) 

+ ^ + ^(4x^ + ^2) 
4 

+ mLo{x^ + y^) 



ai 



aiz 



a2 



aiz + a2 \ z 



+ mLo(x2 + y2) 



2 / ^ V 2 



ai 



The potentials admitting G{ai) depending on several ai allow the existence of different first 
integrals of the form U^G{ai). However, since the system is already superintegrable by including 
a single U"^G{ai), all the other first integrals obtained in this way functionally depend on the 
known ones. 

Because Lq 7^ 0, a necessary condition for the extensibility is the presence of a harmonic term 
in the potential V. 



Extensions of the harmonic oscillators 

As an example, we analyze into details the extensions of the isotropic harmonic oscillator (a par- 
ticular case of El, E3, E7, E8) 

Vi = a3{x'^ + y'^) , 
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and of the anisotropic one (a particular case of E2) 

Va = a3{4.x^ + y^). (10) 

In both cases we assume as 7^ to avoid a trivial potential. 
The only possible extension for Vi is 

= ^ (p' +«3 + / + , (11) 

with the constraint 03 — mLg = due to (3) that sets the potential in the tabulated form. 
The Hamiltonian (11) represents an anisotropic oscillator in and shows explicitly its super- 
integrability being m an integer. Recalling that U = pu — m'^uXi, and by setting X^G = 
a\Px + a2Py = P, we obtain as an example the expression of U'^G 

U'G = Gpi - uplP - \a,GuVu + ^uhuP + f^Gu\ 

The anisotropic oscillator (10) admits, instead, two extensions, corresponding to the two 
different functions G and two different relations between 03 and Lq (items {iv) and {v) of the 
table with ai = 02 = 0). In order to analyze this case in full generality, an iterated procedure 
of extension can be applied. By starting from the one-dimensional oscillator 

Hi = ^pf+ooxj, 
we build a first extension 

H2 = -P2 +Hi^ 2 ^2 = o (^'l + PV + ^^l 2^2' 

z m| 2 mf 

where we use X2 = u and, because of (3), we have Gi = aixi and Lq = uj/mi. The resulting 
potential is an anisotropic oscillator and our procedure produces the third first integral of 
degree mi in the momenta U^^Gi = {p2 — m'[^XHi)"^^Gi. The potential of H2 coincides 
with Va for 772-1 = 2 and to = 4q3 (this is the unique case with a third quadratic first integral). 
A further step is the research of an extension of H2 

Tj 1/2, 2, 2\, 2, ^2, -^0 2 

H?, = - [Pi+P2+Pz) +^^1 H 2^2 H 2:3, 

z m| m2 

where the value of Lq has to be determined. In this case, since the general G2 is (8), condition (3) 
becomes 

aiujxi + 02-^x2 = m2Lo(ao + aixi + 02X2). 
Hence, we get ao = 0, and 

ai ioj - m2Lo) = 0, ^2 ( - "2-2-^^0 ) = 0. 



If ai = 0, we need 02 7^ and 

Lq = — 2 ' ^2 = 023^2- 

m\m2 
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If ai 7^ and 02 = 0, then we have 
Lq = — , G2 = aixi. 

FinaUy, if both ai 7^ and 02 7^ the two conditions are satisfied iff ml = 1 (i.e., H2 is the 
Hamiltonian of an isotropic oscihator) and 

Lq = — , G2 = aixi + 02X2. 
m2 

In the first two cases, H3 represents an anisotropic harmonic oscihator, which is superintegrable 
because mi, m2 are integers, but they have different functions G2 and consequently different 
first integrals U"^'^G2 = {ps — ml^^ Xh2)^^G2- If we set mi = 2, m2 = m and oj = 4a3 in order 
to restrict ourselves to the potential (10), we obtain two possible extensions: if ai = we have 
the relation as — mL^ = that gives the first form in the table. If, otherwise, 02 = we have 
the relation 4a3 — mLg = that gives the second one. 

The extension procedure can be iterated indefinitely obtaining at the n-th step an n-dimen- 
sional anisotropic oscillator with a complete set of first integrals {U^^Gi, C/™'^G2, . . . , V^^^^Gn) 
of degree (mi, m2, . . . , mn-i), that, together with the n Hamiltonians {Hi, H2, ■ ■ ■ , Hn) make 
the system superintegrable. We remark that the systems obtained in this way are characterized 
by the fact that the frequencies are all integer multiples of one of them. See [4, 9] for additional 
details on superintegrability of anisotropic oscillators. 

4 Superintegrable extensions of 

It is straightforward to generalize to E" the procedure previously applied to E^. Let us consider 
in E" with Cartesian coordinates the Hamiltonian 

L = ^ipl+P2^ hpl) +V{xi,X2,...,Xn). 

The general solution of (2) and (3) are 

G = ao + aixi + 02X2 H h a„x„ 

V = mLo [(xi + x?)^ H h (x„ + x°)^] + F(aiX2 - 02X1, . . . , aix„ - a^xi), (12) 

n 

with the constraint ao = OiX^, where x? € M or C, F is any regular function of the arguments 

1=1 

and Lq ^ 0. The corresponding extension is 

H =^pI + mA{L + Vq) + mLoA'^iu + uof. 

We remark that, as well as in dimension 2, the presence of a harmonic term in 1/ is a necessary 
condition for the extensibility. 

Extensions of the three-body Calogero and Wolfes systems 

We consider the particular case of n = 3. If 0.2 = 0,3 = ai and F{Xi,X2) in (12) is 

F = k (Xf 2 + + {Xi - X2)-^) , 
with k gM, then 

F=k( 

V(x-y)2 (x-z)2 {y-zyj 
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coincides with the celebrated Calogero potential, which is a well known superintegrable system 
(see for example [3] and references therein). If, with the same choice for the Oj, 

F = k ((Xi + X2)-2 + (2Xi - X2)-2 + (2X2 - , 

then F coincides with the Wolfes potential, a three-body superintegrable interaction whose 
dynamic equivalence with the Calogero potential is discussed in [3]. If Lq ^ 0, then the first 
integrals of the extended Hamiltonian for m = 2, 3 and F in the Calogero form are respectively 

U'^G = Gpl - 2AupuP - AA^LoGu'^, 

U^G = Gpl - SAuplP - ISA^LqGu^Pu + GA^Lqu^P, 

where G = {x^ + + + xi + 0:2 + X3) and P = Xi{G) = P1+P2+PZ is the conserved linear 
momentum. 

5 Superintegrable extensions of 

In [1] the complete solution G of (2) for mc = 1 (the constant curvature of a sphere of radius 1) 
is computed in standard spherical coordinates as 

G = oo cos 9 + (ai sin ^ + 02 cos cj)) sin 9. (13) 

Equation (3) becomes 

deVdeG + -^d^Vd^G = 2VG. 
sm^ 9 

Therefore, c = ^ 7^ and, by Theorem 1, the extension of L is 

where we assume without restrictions the constants uq, Lq and Wq all zero, and where /t G M. 
Since \k\ can be multiplied by a positive constant simply by rescaling we can assume, if k / 0, 
|k| = Tin? so that the extensions become 

< = 5fJ + ^ (5 (p? + ^i^pj) + 1') . 

For K E C, the explicit form of the extension follows from Remark 2. 

We remark that the numerical factor m? can be absorbed into L by a rescaling of the coordi- 
nates {9,(j)). Whenever V depends on {9,(j)) only through trigonometric functions, the rescaling 
enlights the existence of discrete (polyhedral) symmetries of L on of order depending on m. 

The superintegrable potentials on with first integrals all quadratic in the momenta are 
determined in [5], where the sphere is intended, as previously, as a complex manifold. The 
nine different superintegrable potentials, up to symmetries in 0(3, C) including reflections, 
are, in Cartesian three-dimensional coordinates {x,y,z) with x = sin ^ cos i;^, y = sm9sm(j), 
z = cos 9 [5] , 



K > 0, 
AC = 0, 
AC < 0. 
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«! 0.2 Oi'iW 
9" ~^ —9 ~^ - Q ? 

ai 

-9 5 

\/a;2 + y2' 

_j_ -,2 z2^y2 _|_ ^2 ^2 

yjy"^ + yjw{z - iy) w{z + iy)' 

Ol ^ 02 _^ 03 
x2 y2 ^2 ' 

where it; = x + iy, -u) = a; — iy, oi, 02, 03 G C. We remark that in (x, y, z) we have 

G = a2X + oiy + aoz, 
with ao, ai, 02 G C. 

By putting the previous expressions for V and G given by (13) into (3), we get a hnear 
homogeneous function in (oj) whose components are hnear homogeneous in (01,02,03). The 
non-trivial solutions (with V and G both not constant) of the equation are 





V 


G 


particular cases of 


i 




Q2 031;; 

y)2 yj3 


aoz 


S2 


a 


Ol 
U)2 


j 02 

ld\/ 3;2 + ?/2 


aoz 


S4 


in 




Ol 
lt)2 


GqZ 


SI, S2, 54, 55 


iv 


Ol 


1 "2^ 

z^\/y^ + -2^ 


a2X 


57 


V 




Ol 


a2X + aiy 


52, 53, 57, 59 


vi 




Ol ^ 02 
x2 y2 


aoz 


59 



The only cases without superintegrable extensions for any combination of parameters are 56 
and 58, that are strictly related. The extensible cases of 59 are all equivalent to (v) or (vi) up 
to permutation of the coordinates. Case (v) can be considered a subcase of (iv), but it is listed 
apart because the corresponding expression of G is different. 



52 V = 

53 V = 

54 V = 

55 V = 

56 V = 

57 V = 

58 V = 

59 V = 



Extensions of S9 

As an example of the extension procedure, we develop the computations of extensions and 
first integrals for the case 59. For the subcase (v) we have V = q , with G = (oisin^ + 
02 coscj)) sinO. For m = 3 we have for the Hamiltonian the first integral 

ct ctn'^ Uj 

U^G = B sin ept - 275ctan^'u cos Bpl - 27C^^ — pi + B cos Octanu pipe 
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^ctanii 9 , „ 9 „„ctan^?/ n „^ctan^ii , 

+ 9C——pIp^ - 27Bctan^usmepupj - 27B——pupl - 27C——pj 

ctan'^ncos^ 2 f ctan^nsin^ ctan^u ctan'^u \ 

~ 275 r^- — pj^pe - bAaslB pu + B —pe + C^- ^p^ , 

sm 6 ^ \ cos 6 cos sm 6 cos^ ) 

where B = a\ sin (/> + 02 cos 0, C = oi cos (j) — a2 sin (t> = 

The second of the subcases of 59 admitting a superintegrable extension is (vi) 



sm^ \ sin^ 6 cos^ , 



with G = cos 6. For m = 2 we have for the Hamiltonian H2 the first integral 

tt2^ r, 2 .sin 6* cos 6* 2 , cos6' r,cos9 

U^G = cos Opi - 4 pBPu - ^^rPe - ^ •) ■ 2 n " ^^T^- 



6 Extensions of TTW-type systems 

Let us consider the Hamiltonian 

1 2 , C A o 



For C = 1, X real and 



^(^2) = 2^^ — + ^^-T — : 

cos^ Ax2 sm^ \x2 



L is a generalization to constant curvature manifolds of the Tremblay-Turbiner-Winternitz 
system (see [7, 8]). We consider the possible extensions of L in dimension three given by 
Theorem 1. Since the sectional curvature of the metric of L is "t-c = x the general 

complete solution G of (2) is 

G = aoC^(xi) + (axSc_{x2) + a2C^(2;2)) S-yJyXx), 
and the extended Hamiltonian has the form 
rr _ 2 , X J 

where we assume for simplicity Lq = = Wq = 0. Equation (3) becomes then 

F' (aiCf (X2) - a2QSc^{x2)) = 2CF {aiS^{x2) + a2C^(x2)) . (14) 

If ai = a2 = 0, then the equation (14) is satisfied for all F, including the TTW potential 
with any value of A. In particular, if A is rational then L is superintegrable together with its 
extensions. 

Otherwise, if ai or 02 are different from zero then the solutions F of the equation (14) can 
be obtained after observing that 

C {aiS(;{x2) + a2C^(x2)) = {aiC({x2) - a2C5'((x2)) • 

Hence, 

1 

F 



(aiC^(x2) - a2CS'^(x2))^ 



12 



CM. Chanu, L. Degiovanni and G. Rastelli 



By differentiating the relation ai5^(x2) + a2C^(x2) = ASi^{x2 + 0^ valid for suitable constants A 
and ^, we have 



a result analogous to the one obtained in [1] for the extension of one-dimensional systems. 
Indeed, when G N \ {0} then L is in the form of an extension of a one-dimensional system. 
This is another example of iterative extension. 

7 Conclusions and future directions 

We have shown how the procedure of extension proposed in [1] can be used to produce new 
superintegrable systems starting from the already known ones, together with their first integrals. 
This procedure allows to extend a number of remarkable systems, including TTW and three- 
particle Calogero systems. Moreover, in some cases the procedure can be performed iteratively, 
thus constructing a family of superintegrable systems in higher dimensions. 

Unfortunately not all the superintegrable systems can be extended through our method, but 
this drawback is balanced by the simplicity and compactness of the algorithm that produces the 
constants of motion. Further studies are in progress to find a more general form of extension 
compatible with a larger number of potentials and to analyze iterative extension in other cases. 
New results about the application to nonconstant curvature manifolds of Theorem 1 have been 
obtained. The problem of applying the extension procedure to quantum systems is not yet 
solved: the first integrals described by Theorem 1 cannot be straightforwardly associated with 
symmetry operators. However, for the quadratic first integrals of type f7^(G) a first quantiza- 
tion procedure has been considered, quite unsuccessfully, in [2] and then a second one has been 
studied, leading in suitable cases to symmetry operators. All these progresses will be presented 
in future publications. 
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